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Abstract. We study the density disturbance of a correlated one-dimen- 
sional electron liquid in the presence of a scatterer or a barrier. The 2kp— 
periodic density profile away from the barrier (Friedel oscillation) is com- 
puted for arbitrary electron-electron interaction and arbitrary impurity 
strength. We find that in presence of correlations, the Friedel oscillation 
decays slower than predicted by Fermi liquid theory. In the case of a spin- 
less Luttinger liquid characterized by an interaction constant g < 1, the 
asymptotic decay of the Friedel oscillation is x~ 9 . For a weak scatterer, the 
decay is even slower at small-to-intermediate distances from the impurity, 
with a crossover to the asymptotic x~ 9 law. 



1. INTRODUCTION 

One-dimensional (ID) interacting fermions have attracted a great deal of 
interest in the recent past spurred by the observation that the bosoniza- 
tion technique developed a long time ago [l]-[4] can be used to study the 
interplay between Coulomb interactions and disorder [5]— [12]. Theoretical 
investigations of quantum wires are also stimulated by the possibility of 
performing experiments probing such systems [13]. The last few years have 
seen a tremendous amount of activity, mostly devoted to transport quanti- 
ties like the conductance [5]— [15]. In this work, we focus on the equilibrium 
electron density distribution of an interacting system with broken transla- 
tional invariance [16]. In the presence of an impurity, conduction electrons 
will rearrange in order to screen the impurity charge. For dilute impurity 
concentration, it makes sense to first study the effect of a single impu- 
rity. Depending on the impurity strength, a barrier of arbitrary transmit- 
tance between zero and one hinders the electron flow. The corresponding 
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crossover between an insulator and a metal can significantly be affected by 
the Coulomb interaction between electrons. This issue touches upon im- 
portant physical effects such as the pinning of a Wigner crystal [17] or of 
charge density waves [18], the breakdown of charging effects with increas- 
ing tunnel conductance [19], or quantum transport in ID heterostructure 
channels [20, 21]. 

In the noninteracting case, the impurity is known to lead to a 2kp- 
periodic decaying disturbance of the ground state density, the so-called 
Friedel oscillation [22]. In d dimensions, it has far away from the barrier 
the asymptotic form 

cos{2k F x + r) F ) 
S P( X ) ~d > (!) 

where r\ F is a phaseshift. In dimensions d > 1, the Coulomb interactions can 
presumably be incorporated by Fermi liquid theory such that the asymp- 
totic x~ d behavior is still valid. Since Fermi liquid theory breaks down in 
ID, the asymptotic decay law of the Friedel oscillation can now be modified 
by interactions. Indeed, as will be shown below, due to reduced screening, 
the asymptotic decay is always slower than the Fermi liquid 1 jx law. 

In order to properly describe Coulomb interactions, one has to specify 
the setup under consideration. If one deals with a ID channel in gated het- 
erostructures (quantum wire) [20, 21], the interactions are usually screened 
due to the presence of metallic gates near the channel. A short-ranged 
interaction leads to the Luttinger liquid model [4] described by a single 
dimensionless interaction constant g. The noninteracting Fermi liquid case 
corresponds to g = 1, and the presence of (repulsive) Coulomb interactions 
implies g < 1. On the other hand, one might as well consider a clean iso- 
lated channel where the 1/r tail of the Coulomb potential is not screened 
[23]. We will refer to this situation as the "long-ranged case" in the fol- 
lowing. For the sake of clarity, we shall only consider a single transport 
channel, and most of our analysis will be concerned with the simplest case 
of spinless fermions. Electron-electron backscattering effects can then be 
incorporated by a renormalization of the interaction parameter g [24]. 

Qualitatively, Coulomb interactions favor a smooth density profile. Thus 
one expects that the Friedel oscillation should be smoothed and decay more 
slowly than predicted by Eq. (1). Furthermore, since the Coulomb inter- 
actions lead to a vanishing ground-state conductance in the presence of 
an impurity [5], one should expect a more efficient pinning of the Friedel 
oscillation, i.e., a larger amplitude than in the noninteracting case. Our 
results corroborate this simple picture. For the case of a spinless Luttinger 
liquid, the 1/x law is changed to an asymptotic x~ 9 decay far away from 
the barrier [16]. For a weak scatterer, the behavior of the Friedel oscilla- 
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tion is quite complicated. The decay is even slower than x~ 9 for small-to- 
intermediate distance from the barrier, with a crossover to the asymptotic 
x~ 9 law at some scale xq ~ X^ 1 ^ 1 ^ 9 ^, where A denotes a dimensionless 
impurity strength introduced below. These results have intrinsically many- 
body character and cannot be found from Hartree-Fock or related tech- 
niques. 

In the presence of two impurities new features arise due to the possi- 
bility of resonant tunneling through the double barrier structure [25]. For 
the noninteracting case, it is easily shown that the Friedel oscillation on 
resonance decays faster than x^ 1 , namely with a x~ 2 law resulting from 
the interference of both contributions. Interactions will again modify the 
asymptotic properties of the on-resonance Friedel oscillation. 

1.1. MODEL HAMILTONIAN 

We treat the ID interacting electron liquid in the framework of standard 
bosonization [1]— [4]. This approach is appropriate for low temperatures, 
where only excitations near the Fermi surface are relevant. In the following, 
we will mainly discuss the spinless case in detail. The creation operator 
i/j'(x) for spinless fermions can equivalently be expressed in terms of the 
boson phase fields 9(x) and <f>(x), which fulfill the algebra (we put H = 1) 

[9(x),9(x')U = [<Kx),<f>(x')]-=0 

[<i>(x),0(a/)]- = -(i/2)sgn(x-x') . (2) 

Therefore Il(x) = d x (j){x) is the canonically conjugate momentum to 9(x). 
With ^>t(x) = ip+(x) + ipL(x), where the right- and left-moving parts are 
given by 

V4(a0 = \Itt^- exp[±ik F x + i^(<j>(x) ± 9{x))\ , (3) 
V 2ttvf 

the density operator p(x) = ip jf (x)ip(x) is readily found in the form 

p(x) = — + -l=d x 0{x) + — cos[2/c F x + 2y/^9(x)] . (4) 

7T y/TT IT 

This boson representation of the electron density operator is of essential im- 
portance for our work. The three terms follows. (1) The background 
charge density is kp/w. (2) The density due to right- and left-movers re- 
sults in the second term. (3) The last term originates from the mixed terms, 
i.e., from the interference between right- and left -movers. This 2/jp-term 
is responsible for the Friedel oscillation, since it has a nonzero expecta- 
tion value if translational invariance is broken. The bandwidth cutoff uj c is 
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defined as 

LO c = VFkF , (5) 

which is equal to the Fermi energy for the Tomonaga dispersion relation. 

The effective low-energy theory for a clean noninteracting fermion liquid 
is [4] 

H = ^Jdx [u 2 + (d x e) 2 }. 

Neglecting electron-electron backward-scattering processes, the interaction 
among electrons is then described by the density-density interaction term 

H c = — [ dxdx' d x 6(x)U(x - x')d x >e{x') , 

where U(x—x') is the (screened) Coulomb interaction potential. Most of our 
analysis is concerned with the Luttinger liquid case, where one has short- 
ranged interactions. We are then led to the Luttinger liquid Hamiltonian 



H L = Y J dx 
where the interaction parameter 



n 2 + ^{d x ef 
g 2 



g = . = < 1 

VI + Uo/TTVF 

is related to the forward scattering amplitude Uq , such that g = 1 represents 
the noninteracting case. 

Let us now consider an elastic potential scatterer at x = described by 
a potential V{x). It leads to a contribution Hj = J dxV(x)p(x), and for a 
single 5-scatterer, V(x) = nVkp 1 5(x), one finds from Eq. (4) the generic 
form 

Hi = ^¥-d x 9(0) + V cos[2v^F0(O)] . (6) 

The impurity strength V will often be given as dimensionless quantity 

A = nV/iu c . (7) 

Tuning A from zero to infinity corresponds to changing the transmittance 
of the barrier from unity down to zero. The Hamiltonian is then H = 
Ho + He + Hi, i.e., 

H = — [ dx [n 2 + (d x 6) 2 } + — [ dxdx' d x 6{x)U(x - x')d x ,e(x') 
+ ^-^d x 9(0) + Vcos[2^9(0)} . (8) 

KF 



FRIEDEL OSCILLATIONS IN LUTTINGER LIQUIDS 5 



In the Luttinger liquid case the first two terms are replaced by Hp. The 
model (8) has been the subject of many studies in the past few years, 
primarily with regard to conductance computations [5]-[9]. 



1.2. GENERATING FUNCTIONAL 

As is apparent from Eq. (4), one can compute (p(x = y)) from the gener- 
ating functional Z(y) = (exp[2y/Tr i/j,6(y)]). We can first gauge away the 
forward-scattering term ~ d x 6(0) in Hi by the unitary transformation 
U = exp[ig |2 A0(O)/v / 7r]- Then Z{y) becomes 

Z(y) = <exp[2v^%)]> e^* 8 ^ , (9) 

where the average has to be carried out using UHU' 1 , which is just Eq. (8) 
without the forward-scattering term. 

We formally solve for Z by introducing a field q(r), which is constrained 

by 

g(r) = 2^9{x = 0,r) . 

This constraint is enforced by a Lagrange multiplier field A(r), such that 
one has the effective Euclidean action 



\{d T ef + {d x 9) 2 

V F 



Vp f 

S e [9,A,q] = — J dxdr 

+ — [ dxdx'dr d x 9(x, t)U(x — x')d x /6(x', r) 
2ir J 

+ V J dr cos q(r) 

- 2y/^ind(y, 0) + i J dr A(r) [20F#(O, r) - q(r)] . (10) 

The 9 part of this effective action is Gaussian and can therefore be treated 
exactly by solving the classical Euler-Lagrange equation, 

i d 2 9 t i d 2 e 



v F dr 2 g 2 dx 2 vf 
The solution of Eq. (11) is easily found in Fourier space, 

and similarly for A(r). Then Eq. (11) takes the form 

{uj 2 + u\) 9(k, u) = -2^iv F (A(w) - n e~ iky ) , 
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where we have introduced the plasmon frequency 



Wk = v F \k\yJl + U k /irv F (12) 

with Uk = J dx exp(—ikx)U (x) . In the Luttinger liquid case, this is simply 
v 8 \k\ with the sound velocity v s = v F jg. Defining the boson propagator 
functions 



cos(Zcx) 

^9 e -\gux\/v F ( Luttinger l iquid ) j (14) 



F{x,u>) = v F [°° dk co < kx \ (13) 



one finds the solution 

6{x,T) = ^=j ^ e ibJT [k{u)F{x,u) - nF{x - y,u)]. (15) 

Next we have to determine the action corresponding to the classical solu- 
tion. Inserting Eq. (15) into Eq. (10), we find after some algebra 

Sa[A,q] = |^{[M 2 +AHA(-^)]F(0,u;)-2 / uA(-^ J P(2/^) 

- iA(-(jj)q(u)} + V j dr cos[g + q(r)] - ifiqo , 

where qo is the zero-mode of the auxiliary field. Since S c \ is quadratic in A, 
the Lagrange multipliers are simply found by extremization. The result is 

F(y,u) . q{u) 

A{LO) = /U^r— r + l- 



F(0,u) 2F(0,w) ' 

and inserting this into S c \ gives the generating functional (9) in the form 
of an average over the q field. 

In the end, the generating functional takes the form 



Z(x) = W/( x fV™ 2Asgn ^ (exp 



f du , Fix, ll>) 



\l6) 



where the remaining q average has to be taken with the action 
^ r n f duj q(uj)q(—uj) I' , . 

s ^ = l^^h) +v l dTCOS[q ' ,+q(T)h 

For a Luttinger liquid, we obtain therefore 

/OO A, j I, j I f 
-co ^hr 9 kM|2 + V J dT cos[qo + q{T)] ■ (17) 
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Apart from the kinetic energy term, this action corresponds to the action of 
a Brownian particle (q translates into the position of the particle) moving 
in a cosine potential under the influence of Ohmic dissipation [26]-[29]. The 
damping strength is connected with the Coulomb interaction constant g, 
and the height of the cosine potential is related to the impurity strength A. 
The finite mass of the particle can be associated with the bandwidth to c . 
The envelope function in Eq. (16) 



does not depend on impurity properties at all and involves no q averaging. 
It will turn out that this function governs the asymptotic properties of the 
Friedel oscillation. 

2. FRIEDEL OSCILLATIONS 

In this section, we discuss the ground-state equilibrium Friedel oscillation 
induced by a barrier of arbitrary strength A = ttV/lo c in the presence of 
Coulomb interactions. Most of our analysis will focus on the case of a spin- 
less Luttinger liquid, but we mention some generalizations to the spin-i 
situation or to long-ranged interactions. Since the density profile is sym- 
metric around x = 0, we put x > in the following. 

2.1. GENERAL EXPRESSION 

In the case of a spinless Luttinger liquid characterized by the interaction 
constant g, the boson propagator takes the form given in Eq. (14), namely 



and we have to introduce an ultraviolet cutoff to regularize the integral. 
The appropriate cutoff u c is provided by the Fermi energy [see Eq. (5)], 
and employing an exponential cutoff function exp[—u/u c ], we obtain 



with the microscopic lengthscale a = vp/2guj c . Defining a "lattice spacing" 
a = irvp/uj c , one has a = a/2irg. 




(18) 



F(x,u) = | — 7 exp[— |cj|x/f s ] . 



The envelope function (18) is therefore given by 




W{x) = (l + x/a)' 9 , 
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Let us now compute the deviation in the electron density profile, Sp(x) = 
(p(x)) — hp/ft, caused by the presence of the impurity. From Eq. (16), one 
can verify that away from the impurity the slow (k ~ 0) component in 
Eq. (4) is not affected by a potential scatterer. The only space-dependent 
density profile response to the impurity is the Friedel oscillation, 

5p(x)/p = -(1 + x/a)~ 9 cos{2k F x - g 2 X) P(x) , (19) 

where po = kpj'K. There is a renormalization of the noninteracting phase 
shift 7]p = A, which becomes g 2 X in the case of a Luttinger liquid. 

The pinning function P(x) includes the nontrivial q average and takes 
the form 

(20) 

A useful quantity is the pinning amplitude 

P = -5p(x = 0)/po = P(0) , 

and it is also convenient to define the quantity 

Poo = P(x -> oo) = -(cosg )<? • 

From the definition of the pinning function, it is clear that < P(x) < 1 
must be always fulfilled. Furthermore, the pinning function increases mono- 
tonically from Po to Poo- 

For transmittance one (A = 0), the "charge" q is free and P(x) = 
0. For zero transmittance (A — ► oo), the potential Vcosq locks q at odd 
multiples of tt, and P takes its maximal value, P(x) = 1, for all x. In 
that case, we obtain readily that the Friedel oscillation decays as x~ g . 
This result for zero transmittance can also be obtained by open boundary 
bosonization [11]. The slow algebraic decay of the Friedel oscillation cannot 
be reproduced by Hartree-Fock type calculations and is a true many-body 
effect. For instance, this can explicitly be seen by considering the limit 
A — ► oo, where the Hartree-Fock type procedure devised by Matveev et 
al. [14] would predict a x' 1 decay. In the following, we shall discuss the 
properties of the pinning function in some detail. 

2.2. NONINTERACTING CASE: FERMIONIZATION 

Let us first discuss the exact solution of the bosonized model (8) for the 
noninteracting case g = 1. For this particular value, we can obtain the 
exact solution for all quantities of interest by means of fermionization. This 
is seen by re-writing the Hamiltonian (8) for g = 1 in terms of the right- 
and left-moving fermion operators ip±(x) as given in Eq. (3). The impurity 



P(x) = — ( cos 



Qo + 



f CLLO 



-g\u>\x/v F 
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term is a product of fermion operators at x = 0, and the bulk term becomes 
the massless Dirac Hamiltonian. In the end, H is equivalently expressed in 
the fermionized form 

/d 
dx P^ P {x)-^%{x) + v F \ ^J(o)Vy (0) , 

p=± ' ' ' p,p'=± 

where A = ttV/uj c is the impurity strength (7). 

Since is quadratic in the fermion operators, it is sufficient to study 
the equations of motion 

{^¥t T it) = iA *(*)(40M) +^-(m)) • (2i) 

Away from x = 0, the solutions of Eq. (21) are simply free waves ~ 
exp^—iv^kt ± ikx), where k > is the wavevector. At x = 0, the right- 
and left-moving components are discontinuous. Adding the two equations 
(21) and integrating over an infinitesimal region around x = 0, we obtain 
the jump condition 

01 - (o + ) - o; - vi) (o-) = -2z\ (vi + vi) (o) , 

which implies that the ip± are given in terms of the usual scattering waves 
[30] with /c-independent transmission and reflection amplitudes, or equiva- 
lently a /c-independent phase shift 

Vk = Vf = A ■ (22) 

In the absence of correlations, all effects of the impurity are contained in 
this phase shift, and the Priedel oscillation is readily evaluated in closed 
form. By expressing the reflection amplitude r& in terms of the phase shift 
r/fc, and then using the relation [30] 

Sp( x ) = — dkr k e 2ik \ x \ , 
vr Jo 

we obtain the ground-state result 

Sp(x) = [cos(2A: F |x| + r/ F ) - cos(r/ F )] . (23) 

2tt\x\ 

As expected for the noninteracting case, the ID Friedel oscillation indeed 
decays as 1/x. 
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2.3. SELF-CONSISTENT HARMONIC APPROXIMATION (SCHA) 

Let us begin our discussion of the pinning function for g < 1 by describing 
a simple approximation based on Feynman's variational principle (self- 
consistent harmonic approximation, SCHA) [28, 31]. The most important 
approximation made in the SCHA is the neglect of tunneling transitions 
between different wells of the impurity cosine potential which seems rea- 
sonable for large A. Therefore, we can assume qo to be an odd multiple of 
7r and consider a Gaussian trial action 



StAq] 



f°° dco \u\ . . .. 2 fir 2 , , 



(24) 



where the frequency fi is determined from a variational principle for the 
free energy [32]. It states that the free energy F obeys the inequality 



F < F tT + (H — H tT ) 



tr 



where the average has to be carried out using the trial action (24). Mini- 
mization of F yields 

n = Vexp\-(q 2 ) tI /2 . 



From Eq. (24), we can read off 



(q{u)q{-uj')} tI 



2tt5(uj' - to) 

fi+ M/27T0 ' 



(25) 



such that 



(q 2 ) tr = 2<?ln (l + 
and a self-consistency relation follows, 



2-n-gn J 



Q/V 



1 + 



2itgQ.) 



(2 ff A)f/( 1 -f) 



A > 1 
A< 1 



(26) 



For a strong scatterer A = ttV/u) c S> 1, the trial frequency is simply V, as 
follows by a direct expansion of the impurity cosine term around the minima 
qo, which are odd multiples of n. In this limit, only small fluctuations about 
these minima are possible, with interwell transitions being forbidden by an 
exponentially small WKB tunneling factor. Therefore, we expect SCHA to 
be most valuable for the strong-scattering limit. 

From Eq. (20) together with Eq. (25), one can easily evaluate the now 
Gaussian average. SCHA yields for the pinning function, 



P(x) = exp \-g e ( x + a )/ x o Ex{{x + a)/x ) 



(27) 
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with the exponential integral E\{y) [33] and the crossover scale xq given 
by 

/ 1 /t^/on f V(2sA) , A>1 , no . 

xo/« = ^ (v/n) = | ^V-,) ; A K< x (28) 

where we have used Eq. (26) in the second step. 

Since in the strong-scattering limit xo is even smaller than a, the term 
"crossover" is not meaningful in this limit. Using asymptotic properties of 
E\{y), Eq. (27) becomes for x » max(a,xo) 

p = e -gx /x _ 1 _ yff) 

In the strong-scattering limit, the pinning function is essentially unity for 
all x, and the x~ 9 decay is always found. This result is in accordance with 
Monte Carlo results. 

In the weak-scattering limit, A <C 1, the pinning function exhibits far 
more structure. The crossover scale goes to infinity as A — > 0, namely 
xo ~ A _1 ^ 1 ~ 9 - ) . One expects that there are two different types of behav- 
ior for x <C xo and x » xo- Unfortunately, as we will see below, SCHA is 
unable to provide correct quantitative results except for very strong interac- 
tions, g <C 1. That SCHA becomes more accurate for stronger interactions 
can be rationalized as follows. The presence of interactions leads, loosely 
speaking, to a renormalization of the barrier height. Using a perturbative 
renormalization group (RG) approach [5], one finds that A grows under 
the RG transformation. For strong interactions, it flows quickly into the 
strong-scattering limit, where SCHA is essentially exact. 

For x » xo, SCHA always gives P ~ 1 according to Eq. (29). This is 
an incorrect result, as can be seen from the exact result for the special case 
3 = 1. That failure is due to the complete neglect of interwell tunneling in 
the SCHA. Without tunneling transitions, qo must be an odd multiple of it, 
and one finds Pqo = 1 as predicted by SCHA. However, taking into account 
excursions to neighboring wells, it is readily seen that in general -Pqo < 1. 

Despite of these shortcomings, the effective Gaussian treatment indi- 
cates that for a weak scatterer, there is a crossover, with a slower decay of 
the Friedel oscillation at small-to-intermediate distances x <C xo than the 
asymptotic x~ 9 decay. In fact, Eq. (27) gives 

P(x) = [(x + a)/x ] 9 , x <C x , 

which would imply that the Friedel oscillation does not decay at all up to 
the scale xo- Only for x S> xq, one would have the x~ 9 law. 

In conclusion, for a weak scatterer, SCHA yields for the Friedel oscilla- 
tion 

6p(x)/ P0 = - cos(2fc F x - g A) x j % >> ^ (30) 
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This prediction is certainly incorrect for weak Coulomb interactions, see 
Eq. (23). The effective Gaussian treatment is only valuable for a strong 
scatterer or for strong Coulomb interactions. Furthermore, SCHA provides 
an estimate for the important crossover scale Xq. 



2.4. PERTURBATION SERIES FOR THE PINNING FUNCTION 

To investigate the weak-scattering limit, we next attempt to evaluate the 
pinning function by perturbation theory in A. The perturbation series is 
found by expanding the impurity propagator [27], 



exp 



V fdT cosg(r)| = jr(-V/2) m f P m r^exp 

J J m=0 J {a} 



where one has to sum over all auxiliary variables tij — il? and J 
denotes a time-ordered integration over the m possible intermediate times 
Tj. Thereby, the q average becomes Gaussian again, and P(x) takes indeed 
the form of a power series in A. One finds easily that only odd powers in A 
contribute to the perturbational expansion. 
We find for the lowest-order contribution 



P(x) = \ 7 V(x/a) 1 -° + 0(\ 3 ), (31) 
with the numerical prefactor 



t« = ;:, 2y , (32) 



(1) = 4*-rG7-I 

where T(z) is the Gamma function. This perturbative result is only valid for 

g > \, since otherwise jjp is not defined. That the computation of P(x) is 
indeed a nonperturbative problem for all g < 1 can be seen by computing 
the higher-order terms in A. From dimensional scaling, the perturbation 
series must have the form 



P(x) = ]T A 2 — ) (x/a 



,(2m-l)(l- 9 ) 



Since the higher-order terms increase faster, the first-order estimate (31) 
can only be valid for x <C xq. The crossover scale xo may be computed by 
equating the m = and m = 1 components, or by arguing that P{x) < 1. 
Apart from numerical prefactors of order unity, these estimates for the 
crossover scale coincide with the SCHA crossover scale (28). 
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In contrast to Eq. (30), perturbation theory would therefore predict for 
the Friedel oscillation in the case of a weak scatterer A <C 1, 

5 P (x)/ Po = - C os(2k F x-g*\)J X fW, a l~ 29 ' °<*«*o (33) 

Here, we have assumed that the pinning function approaches a constant 
value Poo < 1 for x 3> xq. This will be confirmed below by Monte Carlo 
simulations. The perturbative result (33) is expected to hold at least for 
weak interactions, 1 — g <C 1. Perturbation theory breaks down even for 
small distances, x <C xq, in the case of strong interactions, g < \. 

In conclusion, perturbation theory predicts a similar weak-scattering 
scenario as SCHA. The Friedel oscillation exhibits a slower decay ~ x l ~ 2g 
at small-to-intermediate distances from the barrier, with a crossover to 
the asymptotic x~ 9 law. The estimate for the crossover scale xq coincides 
in both approximations up to a numerical constant of order unity. From the 
scaling xq ~ \~ 1 /0--a) ^ one observes that the limit of a vanishing barrier 
implies a nontrivial long-distance behavior of the pinning function. 



2.5. MONTE CARLO RESULTS 

One can compute the pinning function (20) for any barrier height A and 
arbitrary interaction constant g by employing numerically exact quantum 
Monte Carlo simulations. As our results for the strong-scattering limit cor- 
roborate the SCHA prediction P(x) ~ 1, we only present numerical data 
for the weak-scattering case here. Numerical simulations are particularly 
useful for a determination of the behavior of the pinning function at small- 
to-intermediate distances from the impurity, x <C xq, where we expect to 
find a power law 

P(x) ~ x 6 * . (34) 

The two approaches discussed above suggest a A-independent exponent 
5 g = g (SCHA) and 5 g = 1 — g (perturbation theory). The SCHA estimate 
should be valid for strong interactions, g <C 1, while perturbation theory 
should hold at least for weak interactions, 1 — g <C 1. 

To compute the pinning function by Monte Carlo, one has to consider 
finite temperatures, such that the frequency integrals become sums over 
Matsubara frequencies, u n = 2irn//3, with (3 = l/fc^T. We have checked 
that the temperatures used in our simulations were low enough to ensure 
that one is in the zero-temperature limit. Typically, (3uj c = 1000 was suf- 
ficient. We have employed a hard cutoff scheme by keeping only the Mat- 
subara frequencies with \co n \ < uj c . The Matsubara components q(to n ) are 
then sampled according to the action (17) using the standard Metropo- 
lis algorithm. Our data were obtained on an IBM RISC 6000/Model 590 
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Figure 1. Monte Carlo data for the pinning function P(y) at g = ^ f° r t\v° different 
values of the impurity strength. The dimensionless space variable is y = uj c x/vf- Notice 
the logarithmic scales. 




Figure 2. Numerical estimates for the exponent S g . Error bars were obtained from 
comparing data for different A. The dashed curves are the SCHA estimate 8 g = g and 
the perturbative estimate 5 g = 1 — g, respectively. 
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Figure 3. Monte Carlo data for the pinning amplitude Po at g = \ as a function of the 
impurity strength A. The dotted curve is a guide to the eye only. 



workstation using 50.000 samples for a given parameter set, with subse- 
quent samples separated by 5 passes. One can compute the full pinning 
function for all x of interest in a single simulation run. 

In Fig. 1, we show data for g = ^ and the two impurity strengths 
A = 0.1 and A = 0.2. The data clearly display the crossover. We find the 
power law (34) for x <C xq, with an exponent 5i/ 2 = 0.28±0.03. For x » xo, 
the pinning function is essentially constant. The asymptotic decay of the 
Friedel oscillation at large x » xq is therefore x~ 9 , while the decay is slower 
for x <C xq, 



6p(x) 



~ x 



-g+Sg 



As is apparent from Fig. 1, the region where this law is valid shrinks rapidly 
as the impurity strength is increased. Our data are generally consistent with 
the scaling xq ~ X^ 1 ^ 1 ^ 9 ^. 

Fig. 2 displays numerical values for S g . From our data, the power-law 
exponent 5 g appears to be independent of A. The SCHA prediction 8 g = g is 
valid asymptotically for g — > 0, while for 1 — g <C 1, the perturbative result 
S g = 1 — g is reproduced, with a smooth turnover between both limits. 

The ability of the scatterer to pin the charge density wave can conve- 
niently be measured in terms of the pinning amplitude Po- As expected, 
Po increases with increasing impurity strength A, see Fig. 3. To study the 
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Figure 4- Monte Carlo data for the pinning amplitude as a function of the interaction 
constant g at fixed impurity strength A = 7r/10. The dotted curve is a guide to the eye 
only. 



dependence of P$ on the interaction strength g, we show in Fig. 4 data for 
Pq at fixed impurity strength. Clearly, the pinning is more efficient as the 
interaction becomes stronger. As mentioned in the introduction, this can 
qualitatively be understood in terms of the scaling of the barrier height [5] . 
Finally, in Fig. 5, we show the quantity Poo for g = \ as a function of A. 
It shows the same qualitative behavior as in the noninteracting case. Most 
importantly, we generally find Poo < 1, in contrast to the SCHA predic- 
tion. Deviations from unity directly reflect tunneling transitions between 
different wells of the cosine potential. 

One might wonder about computing the pinning function exactly at 
the special value g = \- Unfortunately, the fermionization method which is 
able to yield exact results for the conductance [26] cannot achieve this aim 
(except for the quantity Pq), since the density operator (4) cannot be ex- 
pressed in terms of the fermions employed to diagonalize the g = \ model. 
In contrast, for vanishing interaction at the point g = 1, the fermioniza- 
tion method allows for an exact computation of the pinning function. To 
substantiate our results it would be of much interest to apply perturbed 
conformal field theory methods together with the thermodynamic Bethe 
ansatz to the interacting problem. 



0.75 



f£ 0.5 



0.25 



FRIEDEL OSCILLATIONS IN LUTTINGER LIQUIDS 17 



0.75 



0.5 



0.25 







1 


""♦ 


1 1 














i ■ 




g= 1/2 



0.1 0.2 0.3 0.4 0.5 

A 



Figure 5. Monte Carlo data for at g — \ as a function of the impurity strength. 
The dotted curve is a guide to the eye only. 



2.6. LONG-RANGED INTERACTION 

In the absence of metallic gates, one has to take into account the long- 
ranged character of the Coulomb interaction between the electrons in the 
ID channel. The 1/x tail of the potential leads to a k — > divergence of 
the Fourier transform, 



00 cos(kx) 
dx- 



ly/d 2 + 



(2e 2 / K ) K (kd) ~ (2e 2 /K) \ In | kd\ \ 



«1, 



where d is the width of the ID channel (kpd <C 1), and k denotes the 
dielectric constant. The Bessel function K (z) can be approximated by a 
logarithm here, and the corresponding logarithmic correction in the plas- 
mon dispersion relation (12) is 



= a c v F \k\^J\\n\kd\\ , 
where the dimensionless Coulomb interaction constant 



a c = J 2e 2 /ttkvf 
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Figure 6. Monte Carlo data for Po as a function of impurity strength in the long-ranged 
case. We have taken a c = 2.5 and d = -kvf/lo c . 



is of order unity in typical quantum wires [21]. 

Unfortunately, the boson propagator (13) cannot be evaluated in closed 
form anymore. However, one can obtain the asymptotic expansion for large 



F(x,u) 



2v F 



■ sm 



\UJ\X 



(35) 



_a c v FV / \ln\uj\ 

Remarkably, up to a prefactor of order unity, this expression holds even at 
x = 0, such that [9] 

F(0,u) 



1 



to 



UJ 



Computing the pinning function by Monte Carlo, we find that P(x) is 
essentially constant, P(x) ~ Pq, as could have been anticipated from the 
corresponding limit g — > of the Luttinger liquid case. The functional 
dependence of Po on the impurity strength is depicted in Fig. 6. 

For long-ranged interactions, the decay of the Friedel oscillation away 
from the barrier is completely governed by the envelope function W(x) 
specified in Eq. (18). Asymptotic evaluation of the frequency integral using 
(35) shows that for x S> d, the leading contribution is 



W(x) ~ exp 



—const 



FRIEDEL OSCILLATIONS IN LUTTINGER LIQUIDS 19 



This decay is slower than any power law. In conclusion, in the presence of 
long-ranged unscreened Coulomb interactions, the Friedel oscillation takes 
the following form far away from the barrier 

5p(x) ~ cos(2/cpx) exp — const. Vinx . (36) 

This long-ranged density disturbance should lead to strong quasi-Bragg 
peaks in x-ray scattering. 

2.7. SPIN-i CASE 

Let us now briefly comment on the spin-i case. The bosonized treatment 
proceeds very similarly [16]. One can introduce charge fields 6 p (x),(pp(x) 
generalizing the fields 6{x) and 4>(x) employed previously, and, in addition, 
spin fields a (x) and (j> a (x). These fields are linear combinations of the re- 
spective spin-up and -down fields, such that the electron creation operator 
for spin s = ± at location x takes the form 



27rv ? P =± 



ip k F x + Jn/2 [9 p {x) + s9 {x)\ 



x exp 



iJn/2 [<f>p(x) + scfrvix)] 



The various 6 U fields {y = p, a) commute among themselves, as do the § v 
fields. The correct generalization of (2) is then given by 

\4> v {x),d vl (x')\- = -{i/2)8 vv i sgn(x - x') , 



J x y rv 



and the canonical momentum for the 9 U field is Ii u = d x 

We are concerned with density distributions in the presence of impurities 
or barriers. Generalizing Eq. (4), the bosonized form of the density operator 
for spin-| electrons is 

p(x) = p + J2/ird x 6Jx) + — cos[2k F x + V2tv9 p (x)] cos[V2tt 6 a (x)] 

V 7T 



+ const. cos[AkFX + v 8ir p (x)] , 

where the background charge is now po = 2kp/TT. The three other terms 
are (1) the long-wavelength contribution, (2) the 2k F charge density wave 
part, and (3) the 4fcp Wigner component [3, 23]. The Wigner component is 
not present in the spinless case, since two right-movers have to be flipped 
into left-movers simultaneously for this term to arise. Due to the Pauli 
principle, this can only happen for spin-^ electrons. 
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Assuming to be away from lattice or spin density wave instabilities, 
and neglecting electron-electron backscattering for the moment, the clean 
system is described by H = H a + H p with 



dx 



U 2 p (x) + (d x 0„(x)) 2 ]+i / <i.nlx' <),(>,< .r)U{.v-.v')i). t ,0 p {.v') 



The spin part H a is identical to the charge part H p with no interaction 
potential U and the p fields replaced by the a fields. Furthermore, the 
impurity backscattering contribution has the form 



Hi = Fcos[V27r<9 p (0)] cos[V2tt 6> ff (0)] . 

Spin and charge parts are now coupled through this term. The forward- 
scattering contribution ~ d x 6 p (0) leads only to a phase shift and is readily 
included by a gauge transformation. 

Friedel oscillations can again be extracted from a generating functional 
generalizing Eq. (9), 



Z(x,p p 



cxp 



2iri ^ Vu°u( x ) 



iy=p,cr 



Since the impurity influences 9 U only at x = 0, we constrain the field 
amplitudes 9 v {x = 0) to be equal to new fields, q v {T) = \/2k 0„(O, r). 
We then proceed as before, and in the end, one is left with the nontrivial 
average over the q v fields alone, which are coupled to each other through 
Hi. 

Collecting together all terms, we obtain 



z=-v{x) n wfi' 2 { 

u=p,a 



x) 



(37) 



The envelope functions W v {x) generalizing (18) are 



W u (x) = exp 



2vr 



FH(0, 



The charge and spin boson propagators are defined as in Eq. (13), and for 
a Luttinger liquid, we have again a single interaction constant g = g p < 1. 
In the absence of a magnetic field and any spin-dependent interactions, one 
has to put g a = 1 to respect the SU(2) spin symmetry [4, 5]. Thus, the 
boson propagator functions read 



F (ly) (x,w) = exp[-\g u ux\/v F ] . 
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Finally, the quantity V in Eq. (37) generalizes the pinning function. We 
find 



I tt ( f dw ( x 

Hx) = -( II ex P ( Q»,o + J ^ q ^JF(^ 

\V=p,<7 \ L ~ ^ ' 



The q bracket stands for an average taking the action 
duo M 

/ 

v=p,a ' 



UJ 



S[<lA = f ■7rJ-^-\qv( UJ )\ 2 + V [ dTcos[q afi + q a (T)}cos[q pfi + q p (T)} . 



In the following, we only discuss asymptotic properties far away from the 
barrier. For these, we can assume that V is constant such that only the 
envelope functions W u (x) have to be evaluated. The spin part is like in 
the noninteracting case and gives rise to a l/\/x factor in the 2kp Friedel 
oscillation. We note that there is the same crossover as discussed above, 
with a slower decay of the Friedel oscillation at x <C xq. 

Combined with the charge channel, we find an asymptotic ~ x ~( 1+9 ^ 2 
decay for the 2/cp Friedel oscillation. This is slightly faster than the corre- 
sponding x~ 9 law for spinless electrons and reflects that one starts to go 
away from the extreme ID case by incorporating the second (spin) channel. 
As electron-electron backscattering cannot be treated by a simple renor- 
malization of g in the spin-| case, there could be weak logarithmic correc- 
tions to the ~ x~( 1+9 ^ 2 law, depending on the magnitude of the backscat- 
tering amplitude [24]. Generalizing the two-channel spin-degenerate case 
discussed in this subsection to a multi-channel situation, one would then 
expect a x~ l law for large channel number, since Fermi liquid theory will 
eventually be valid in that case [7]. 

Remarkably, for spin-i electrons, there is also a 4fcp Friedel oscillation 
component 

Sp(x) ~ cos(4/cpx) x~ 2g , 

which dominates over the 2Jcf contribution for strong enough correlations, 
g < ^. Since 4fcp corresponds to the inter-particle spacing, this suggests 
that for <? < 5 signatures of Wigner crystal behavior are induced by the im- 
purity. The AIzf component has been seen in numerically exact small-chain 
calculations by Hallberg and Balseiro [34]. In their Lanczos calculation, for 
a Hubbard chain with a magnetic impurity, the 4k f component was found 
to be present for strong interactions. 

Wigner crystal behavior has also been found by Schulz [23] for the clean 
system with long-ranged correlations. For 1/x interactions, the 4fcp Friedel 
oscillation decay is extremely slow. While the spin degrees of freedom in- 
volve again the 1/y/x factor suppressing the 2k.F component, the 4kp Friedel 
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oscillations decay again like exp(— cyln~r), i.e., slower than any power law. 
Effectively, one will then only observe the 4kp component. In the spinless 
case, the same quasi long-ranged behavior appears for the 2fcp component 
already because the spin channel is absent, see Eq. (36). 

2.8. FRIEDEL OSCILLATIONS NEAR A DOUBLE BARRIER 

Let us now consider a double barrier arrangement which allows for the 
possibility of perfect transmission even in presence of Coulomb interactions. 
The resonant tunneling problem for a Luttinger liquid has been studied 
extensively before [25], and here we shall focus on the Friedel oscillation for 
a spinless single-channel Luttinger liquid. The symmetric double barrier 
structure considered is described by the Hamiltonian 



instead of the single-barrier expression (6). We omit the forward-scattering 
contribution since it is inessential for the subsequent discussion. The res- 
onance can then be tuned by varying fcpi? (or by varying a gate voltage 
coupling to the charge on the "island" formed between the barriers), and 
we employ again a dimensionless impurity strength A = ttV/lo c . 

To study the Friedel oscillation for the interacting double-barrier prob- 
lem, we proceed as in the single-barrier case and compute the generating 
functional (9). The 9 field can be integrated out by introducing the two 



where these constraints are enforced by Lagrange multiplicator fields A± (r) . 
The Gaussian integration over the 9 degrees of freedom and the ensuing 
minimization of the Lagrange multiplier fields leave us with the following 
expression for the Friedel oscillation 



with po = /cp/7T. The ground-state envelope function takes the general form 



Hi = V cos[2^9(aR/2)] 



(38) 



a=± 



fields 



q± = 2^9(±R/2) , 



Sp(x)/p 



W 2 (x)Ree 2ikF ^P 2 (x) 



W 2 (x) 



expj^{-F(0,u,) + 




x 



F 2 {\x + R/2\,uj) + F 2 (\x - R/2\,u) 
2 F J^ ^ F(\x + R/2\, u)F(\x - R/2\,u) 
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with the boson propagator F(x,u) defined in Eq. (13). For the spinless 
Luttinger liquid, this gives to the right of the right barrier 

W 2 {x) = [l + {x- R/2)/a}~ 9 , x>R/2, 

where again a = v s /2uj c . 

The pinning function P 2 (x) contains all dependency on impurity proper- 
ties and can be expressed as average in {q + , q_} space. It is more convenient 
to switch to the symmetric and antisymmetric combinations [25] 

q = (q+ + q^)/2 
Q = (<7+-<Z-)/2 

describing the transmitted charge q and the island charge Q, respectively. 
The action for the impurity averaging reads for finite temperatures 

1 00 l l 

or ^11 \ UJn \ QnQ-n 



/3 2irg 1 + exp(-|w n |i?/v s 

QnQ—n 



n=—oo 

+ 

1 

n=—oo 



\ - \u} n \ 
3 ^ 2kqT 



(3 n f^ oc 2irg 1 - exp(-\u; n \R/v s ) 

+ 2V / dr cos g(r) cos[k F R + Q(t)\ , (39) 
Jo 



where q n and Q n are Matsubara components at frequency uj n = 2-irn/fl. The 
two modes are coupled by the impurity term. The Q mode has acquired a 
mass gap now [5]. 

For the spinless Luttinger liquid, the complex-valued pinning function 
P 2 (x) takes the general form 

/ r { ( e -\Wn(x+R/2)\/v s + e -\u> n (x-R/2)\/v s 

P 2 (x) = -(exp[-^( i + e -Klfl/«. qn 

e -\w„{x-R/2)\/v a _ e -\u) n {x+R/2)\/v a 
+ 1 _ e -\w n \R/v s Qn )\/q,Q ' 

To the right of the right barrier, x > R/2, this yields the simpler expression 
P 2 (x) = - (exp \W e-\^ x - R ' 2 ^(Q n + q n )] \ 

\ V P n ]/ q ,Q 

where the average has to be taken using the action (39). 

In the following, we shall focus on the case of small barriers. This has two 
reasons. (1) In the strong-scattering limit the resonances are very sharp, 
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and one might not be able to "find" the on-resonance Friedel oscillation. 
(2) A strong-scattering SCHA treatment, as applied to the single-barrier 
case previously, is not possible here since SCHA neglects tunneling and 
hence does not capture resonant tunneling. The SCHA prediction is always 
x~ 9 , which is only correct in the limit A —> oo where the barriers reflect 
completely and no resonance is possible anymore. 

In the weak-scattering case, we may compute the pinning function per- 
turbatively in the impurity strength. To lowest order in A, we obtain the 
Friedel oscillation 



with 73 defined in Eq. (32). This result is only valid for g > ^ and A <C 
1. Off resonance, we find the usual x l ~ 2g dependence, see Eq. (33). On 
resonance, the 2kp component of the Fourier transform = 2V cos(kR/2) 
of the scattering potential (38) vanishes [5]. Hence, the two cosine terms in 
Eq. (40) have different signs and interfere destructively. In that case, the 
Friedel oscillation has a different asymptotic behavior at x S> R, 



8p(x)/p = -A (2g - 1) 7 « cos[2k F (x - R/2)} - - . (41) 

^ rv \ rv I 



The third-order contribution ~ A 3 goes like x~ Ag on resonance. Therefore, 
Eq. (41) holds in the asymptotic regime, contrary to the off-resonance case, 
where the higher order terms decay slower than the lowest-order contribu- 
tion and thus lead to a modification of the perturbative prediction in the 
asymptotic regime (instead of x 1 " 29 , the correct asymptotic off-resonance 
law is x~ 9 ). 

In conclusion, for resonant tunneling through a double barrier, the 
Friedel oscillation can be dramatically affected if the barrier is tuned on 
resonance. The usual x~ 9 law is turned into the faster x~ 2g law far away 
from the scatterers. This should be observable as a much weaker signal in 
x-ray scattering or NMR data. 



In this article, we have applied the bosonization method to compute Friedel 
oscillations for interacting fermions in one spatial dimension. The success 
of this method might come as a surprise, since density properties usually 
depend on details of the (noninteracting) electron band, such as band cur- 
vature. However, the asymptotic properties (which means several lattice 
spacings away from the barrier in that context) depend only on Fermi sur- 
face quantities and can be computed by bosonization. 



Sp(x)/p = -A7W J2 cos ( 2 M^ - [(x - aR/2)/a] 



l-2g 



(40) 




3. CONCLUDING REMARKS 
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We have shown that in ID the Friedel oscillation in a correlated fermionic 
system decays slower than in the Fermi liquid case. This is in conflict with 
naive expectations based on the smeared momentum distribution of a Lut- 
tinger liquid. While in a Fermi liquid thermal smearing of the Fermi distri- 
bution induces an exponential decay of the Friedel oscillation, the smeared 
momentum distribution of a Luttinger liquid does not imply a faster de- 
cay. This shows again the failure of the quasiparticle picture for correlated 
fermions in ID. 

The results provided here have intrinsic many-body character and, to 
the best of our knowledge, cannot be obtained by other methods avail- 
able at the moment. While results for the conductance can be obtained 
for weak interactions from a Hartree-Fock type approach [14], it is not 
possible to obtain the correct density profile with such methods. Further- 
more, the fermionization technique allowing for exact results for the special 
interaction constant g = \ fails for the density profile as well. We hope 
our findings will motivate further theoretical work, e.g., on the effects of 
magnetic impurities on the screening cloud in one-dimensional interacting 
fermions. 
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